Appearance of Light Clusters in Post-bounce Evolution of Core-Collapse Supernovae 
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We explore the abundance of light clusters in core-collapse supernovae at post-bounce stage in 
a quantum statistical approach. Adopting the profile of a supernova core from detailed numerical 
simulations, we study the distribution of light bound clusters up to alpha particles (2 < A < 4) as 
well as heavy nuclei {A > 4) in dense matter at finite temperature. Within the frame of a cluster- 
mean field approach, the abundances of light clusters are evaluated accounting for self-energy, Pauli 
blocking and effects of continuum correlations. We find that deuterons and tritons, in addition to 
''He and ''He, appear abundantly in a wide region from the surface of the proto-neutron star to the 
position of the shock wave. The appearance of light clusters may modify the neutrino emission in 
the cooling region and the neutrino absorption in the heating region, and, thereby, infiuence the 
supernova mechanism. 

PACS numbers: 21.45.-v,21.65.-f,26.50.+x,97.60.Bw 



The formation of clusters in nuclear matter is one of 
pivotal phenomena in nuclear structure and heavy ion 
collisions. Light clusters such as "^He appear also in as- 
trophysical environments at certain densities and tem- 
peratures. If they exist abundantly, they may change 
the properties of dense matter and contribute to new re- 
actions controling the dynamical evolution of stars. 

In the explosion mechanism for core-collapse super- 
novae, which has been a long standing problem despite 
continuous research for decades [l| , the existence of "^He 
and nuclei in the heating region has been proposed to be 
a possible agent to help the stalled shock wave via the 
neutrino heating mechanism 0, . In most of the recent 
numerical simulations of core-collapse supernovae study- 
ing the gravitational collapse of massive stars, the shock 
wave produced by a core bounce stalls on the way and 
needs some mechanism to revive the outward propaga- 
tion of the shock wave. The heating source through neu- 
trino reactions in dense matter behind the stalled shock 
is one of key issues in sophisticated numerical simulations 
It has been shown that neutrino-^He reactions may 
marginally affect the standing accretion shock instability 
and consequently the supernova shock revival fG] . 

Besides ■*He, the light clusters such as deuterons, tri- 
tons and '^He have not been included in numerical sim- 
ulations of supernovae so far. Those light clusters may 
become new agents for absorption and emission of neutri- 
nos, which can carry energy away and heat the material 
in the important phase of supernova dynamics. In the 
currently available sets of equations of state (EOS) for 
supernova simulations such as the Shen-EOS ^7i| and the 
Lattimer-Swesty EOS f^, dense matter is described as a 
mixture of neutrons, protons, *He, and one representa- 
tive species of nuclei. Since the occurrence of light clus- 
ters in addition to *IIe is a natural outcome in nuclear 
statistical equilibrium, one should clarify whether they 
are also present in realistic profiles of a supernova core. 



A consistent description of the composition of hot and 
dense nuclear matter is preferable for the rigorous treat- 
ment of neutrino transport in the post-bounce evolution 
of core-collapse supernovae. Progress in the description 
of clusters in dense matter [9|, I lOl based on a generalized 
Beth-Uhlenbeck approach [II], [I^l enables us to evaluate 
the abundance of deuterons, tritons, and helium nuclei 
altogether within a microscopic approach. The mixture 
of tritons and ^He and their neutrino reactions have re- 
cently been studied under the condition of fixed temper- 
atures and proton fractions ^IS^, but without deuterons 
and heavy nuclei. 

In this paper, we report the appearance of light clusters 
in dense matter of supernova cores adopting the micro- 
scopic treatment of clusters developed within a quantum 
statistical approach. We aim at the unification of the 
quasiparticle approach such as the relativistic mean-field 
theory 0, which may be adequate at high baryon 
densities, and the cluster approach including the virial 
expansion [ill, [13, [Hj [Hj [3 which may be applicable in 
the low-density limit. 

The EOS of nuclear matter ni{T, ^p, /i„) giving the to- 
tal particle densities of protons and neutrons, i = p,n, 
respectively, as function of temperature T and chemi- 
cal potentials /i,;, is obtained from the spectral function 
A(l,w), with 1 = {fci,(Ti,Ti} denoting momentum, spin 
and isospin of a single nucleon. In the low-density, low- 
temperature region, a nonrelativistic Hamiltonian ap- 
proach is possible which is based on a nucleon-nuclcon 
interaction F(12,1'2'). Well known examples are the 
Bonn and Paris potentials or its separable representa- 
tions, reproducing empirical properties such as scattering 
phase shifts. Within this quantum statistical approach, 
the baryon density is calculated as 
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Further thermodynamic quantities such as pressure 
p{T, fj,p, or free energy density /(T, n^, n„), being po- 
tentials, are obtained by integrations. 

The spectral function is related to the self-energy E, 
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Besides the quasiparticle contribution at £^'^"(1) — 
p\/2mi + Rel](l, it contains also the contribu- 

tion of bound states and continuum correlations included 
in ImS](l,Lj — iO), so that within a generalized Beth- 
Uhlcnbeck approach a virial expansion can be performed 
in the low-density region. After the cluster decomposi- 
tion of the self-energy [ll|, [l^ , the contributions of the 
different clusters are treated in a systematic way. 

The cluster decomposition for the self-energy involves 
the calculation of the A-nucleon T matrix which obeys 
a corresponding Bethe-Salpeter equation. Two-particle 
correlations in nuclear matter can be accounted for by 
considering the two-particle Green function in ladder ap- 
proximation. The corresponding Bethe-Salpeter equa- 
tion contains self-energy and Pauli- blocking terms [l^. 
The cluster-mean field approximation ^15j yields a wave 
equation which describes also arbitrary clusters in a clus- 
tered medium. For the two-particle system, neglecting 
correlations in the surroundings, the following wave equa- 
tion results 
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The medium Fermi distribution function /(I) = 
{exp[£;'i"(l)/T - fix/T] + contains the effective 

chemical potential (li which is determined by the total 
proton or neutron density calculated in quasiparticle ap- 
proximation, rii = "Y^^ f{l)S-r-^^i. It describes the occu- 
pation of the phase space neglecting any correlations in 
the medium. For the parameter values T, ns considered 
below the account of correlations in the medium gives no 
significant changes in the spectral function [l5l | . 

The quasiparticle energies E'^^{1) are self-consistently 
calculated if we evaluate the self-energy or the corre- 
sponding T matrix. Using the cluster decomposition 
and restricting only to the two-particle contribution , 
the results coincide with the Brueckner K matrix includ- 
ing hole-hole contributions. The quasiparticle energies 
can be parametrized and directly related to empirical 
data on finite nuclei. An appropriate parametrization 
can be found on the basis of relativistic mean-field the- 
ories 0, [3| which also involves the results of relativistic 
Brueckner-Hartree Fock calculations [l3|- We will use 
the TMl parametrization here. 

From the solution of this in-medium two-particle 
Schrodinger equation ^ the scattering states and possi- 
bly the bound states are obtained. Due to the self-energy 
shifts and the Pauli blocking, the binding energy of the 



deuteron Ed{P; T, ^p, fin) as well as the scattering phase 
shifts 6\{E, P; T, fip, ijl„) in the respective isospin singlet 
or triplet channel A, will depend on the temperature and 
the chemical potentials [l^, [l^ . 

In-medium Schrodinger equations similar to Eq. ([2]) 
containing quasiparticle shifts and Pauli blocking terms 
are derived for higher clusters with mass number A 
and charge Z [l^. The shift of the bound state ener- 
gies £;X^(P;r,nO = Ea^z + Af;^ + A^^' + AS°^i, 
containing in addition to the single-particle self-energy 
shift A^^2 and the Pauli blocking term Aj^**^'' also the 

Coulomb shift A5™\ can be calculated within perturba- 
tion theory. Besides the quasiparticle shift at zero mo- 
mentum A^J^^ = ZE^''{0) + {A-Z)E^''{0) which can be 
included into the chemical potential, there is the contri- 
bution due to the effective mass m*. Assuming for light 
clusters 2 < A < A a Gaussian wave function with the 
nucleonic rms radii {r'^)A,z, the self-energy shift results 
as Af'z = A^'^^° + 3(yl - l)h%\z{^ - m*)/{im^) af- 
ter introducing Jacobian coordinates and separating the 
c.o.m. motion, where h\ z — 3(^4- 1) / {A{r'^) a,z)- With 
Gaussian wave functions we obtain the following estimate 
for the Pauli blocking shift 
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Similar expressions [l9| can be given for the weakly 
bound clusters with 5 < A < 11. The Coulomb shift 
A5°z' is calculated in Wigner-Seitz approximation, see 
also [l^ . Within the parameter values considered be- 
low, the influence of the Coulomb corrections on the com- 
position is small. 

For heavier clusters with mass numbers A > 12, the 
self-energy and Pauli-blocking shifts become less impor- 
tant and will be neglected here. The heavier clusters 
repel the nuclear matter so that the mean-field effects 
are mostly produced by the other nucleons within the 
cluster and are contained in the cluster binding energy. 
For a more detailed consideration see p^ . 

Using the cluster decomposition of the self-energy, 
from Eq. ^ a generalized Beth-Uhlenbeck formula can 
be derived which is formally similar to the virial expan- 
sion of the EOS 



OO 

ns(T,^p,/z„) = ^ 

A=1.,Z 



A nA,ziT,fip,p.n) ■ (4) 



Introducing the Fermi or Bose distribution function 

EOS (|3]) contains the single-particle contribution 
ni^z = '^'12p h,z{E^^{P)), the two-particle contribu- 
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tion of medium-modified deuterons and scattering states 
in the isospin singlet and triplet channel, see [T^.[l8j). as 
well as the contribution of all higher clusters. 



3 



The bound state contribution of the ^-nucleon clus- 
ter to the EOS (H has the form [H [H, [H n^.z = 

Ep,x[fAMEXz,xiP)) - fAAETSiP))] where A de- 
notes the internal degrees of freedom such as spin and ex- 
citation states. The summation over the c.o.m. momen- 
tum P runs only over the region where the in-medium 
binding energy E'^^ (P) is lower than the edge of the 
continuum E'^^'^{P) = A Ei''{P/A). The summation 
over all excited states of a large cluster can be replaced 
by an intergal over the excitation energy using the den- 
sity of states, see [sl. IT9|. 

Of importance is the contribution of continuum corre- 
lations which can be represented by the scattering phase 
shifts. Extended investigations have been performed with 
resp ect to the contribution of two-nucleon phase shifts 
They reduce the contribution of the bound states. 
Partial integration, using the Levinson theorem, gives 
the second contribution fA,ziE'^'f{P)). Thus, there is 
no jump in the EOS when a bound state merges with the 
continuum if, e.g., the baryon density is increased. The 
remaining correlations in the continuum are compensated 
to a large extent by the contribution of interaction to the 
quasiparticle shift and will be neglected. Explicitly this 
can be shown for the contribution of the two-nucleon sin- 
glet and triplet scattering phase shift contributions. 

In order to explore the fraction of nucleons bound in 
light clusters, we utilize a snapshot of a supernova core 
at post bounce stage, which is important for setting the 
revival condition of a shock wave and the resulting ex- 
plosion. We use the profile obtained by the detailed nu- 
merical simulation [20| of the gravitational collapse and 
core bounce for a ISM© pre-supernova star adopting the 
Shen EOS 0. 

Figure [T] displays the profiles of density (p) , electron 
fraction (Ye) and temperature (T) at 150 ms after the 
core bounce. The shock wave is stalled at ^130 km. 
The proto-neutron star is just born at the center and 
matter is falling down onto it through the shock wave. 
The neutrino sphere, where the neutrinos are emitted, is 
located at 20-80 km depending on energies and species. 
The electron fraction is ^0.5 in the iron nuclei region 
at r >100 km and becomes smaller (neutron rich) for 
dissociated nucleons inside due to the balance under beta 
equilibrium including neutrinos. Note that the electron 
fraction is the same as the total proton fraction due to 
charge neutrality. 

At the shock wave, the iron nuclei are dissociated into 
"^He, and, subsequently, nucleons at high temperatures 
due to the shock heating. We remark that the nuclear 
statistical equilibrium is maintained during the dynamics 
because the temperature and density are high at the re- 
gion of current interest and, therefore, the compositional 
change proceeds fast enough via charged particle reac- 
tions. During this process of dissociation, light clusters 
may appear in addition to ^He and can modify the cool- 
ing and heating rates through weak reactions with lep- 
tons. Conventionally, only nucleons are considered for 
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FIG. 1: Density (a), electron fraction (b) and temperature 
(c) profile of the supernova core at 150 ms after core bounce 
as a function of the radius. The position of shock wave and 
the boundaries of heating and cooling regions are indicated 
in the top and bottom panels. 



the cooling process via neutrino emission at the surface 
of proto-neutron star (cooling region) and the heating 
process via neutrino absorption behind the shock wave 
(heating region). 

We solved the EOS dH) for given profiles of T, p, Ye as 
shown in figure [TJ The masses of nuclei are taken from 
Audi and Wapstra |2l|. The chemical potentials p,p,pn 
are determined in a self-consistent way. The abundances 
of the different clusters have been calculated, and frac- 
tions Xa,z = AnA,z/nB of nucleons bound in deuterons, 
tritons, '^He, ^He as well as the fraction Xa of nucleons 
bound in nuclei with A > 5 are shown in figure [51 

The contributions of free protons and free neutrons to 
the baryon density follows nearly the behavior given by 
Shen et al. 0, H^l. In the inner region (r < 30 km) 
where the density and temperature become high, dom- 
inant clusters are deuterons and tritons because of the 
high neutron fraction. Here the main contribution to 
Xa is from A = 5. In the outer region (r > 100 km), 
^He and heavy nuclei become abundant. In the region 
between, the deuteron fraction is relatively high and the 
fraction of ''He rises to about 90 % at r ~ 150 km, be- 
fore it goes down when the heavier clusters, A > 5, are 
formed. The wiggle of Xa near r ~ 150 km is caused 
by the weak binding energy of nuclei with 5 < A < 11 
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FIG. 2; (Color online) Mass fraction Xi of light clusters as 
function of the radius for the post-bounce supernova core 
shown in Fig. 1. 

which are of relevance there. At larger radii, the strongly 
bound nuclei of the iron region will dominate. 

The appearance of light clusters in supernova cores 
at post-bounce stage indicates that there are additional 
weak reactions involving them which will have influence 
on the supernova dynamics. The neutrino interactions 
on light nuclei, which appear inside the proto-neutron 
star at ~10 km, will contribute to the opacity. This 
may lead to a modification of the diffusion coefficient 
for neutrinos, the time scale of neutrino diffusion and 
the location of neutrino sphere, where most of the neu- 
trinos start streaming freely. The light nuclei at ~60 
km in the cooling region are possible sources of neu- 
trino emission through electron and positron captures. 
They will modify the neutrino luminosities and spectra 
because the emission rates per nucleon for light nuclei 
may be effectively smaller than the large emission rates 
for nucleons. Medium effects on the neutrino emissiv- 
ity have been considered in (2^ using the many-particle 
approach described above. In fact, the abundance of nu- 
cleons is reduced by about 10% due to the appearance of 
light clusters and the corresponding neutrino luminosities 
may be lowered. We expect that these modifications of 
the neutrino fiux and energy change the neutrino heating 
rates. 

The light clusters in the heating region behind the 
shock wave become the target compositions for neutrino 
absorption and will modify the heating rate to assist the 
revival of shock wave. The dominance of ^He at ^150 
km in the heating region may contribute to the heating 
rate via neutrino absorption as discussed by Haxton et 
al. 0, Q. It is interesting to see that the fraction of 
^He in the current approach is larger than the value in 
the Shen EOS [13] where the low binding energies for 
the clusters with 5 < A < 11 are not accounted for. 



Deuterons become more abundant than ^He at r < 100 
km and dominate over other nuclear species in the heat- 
ing region. The neutrino absorption cross section on 
deuteron is much larger than that for ^He [2^, [2^ . 
therefore, deuterons may become more effective targets 
for the heating than "^He. Neutrino heating processes 
on triton and '^He should be studied whether they may 
contribute as well [l^. When nucleons appear together 
with light clusters, the heating via the neutrino absorp- 
tion on nucleon is dominant since the absorption cross 
section on nucleon is largest among others. In this case, 
the appearance of light clusters reduces the total heating 
rate. These processes with light clusters will influence the 
neutrino-heating mechanism in the revival of the shock 
wave for a successful explosion. 

We have adopted the spherically symmetric snapshot 
at 150 msec after bounce as an exemplary case. Other 
time sequences and multi-dimensional aspects will be 
studied in our subsequent paper to examine whether 
light clusters appear generally in a significant amount. 
It is known that nuclei appear more abundantly behind 
the shock wave in the hydrodynamical instability and 
the neutrino-^He reactions may contribute to the heat- 
ing 0. It would be interesting to study the fractions of 
light clusters in profiles with wider ranges of density and 
temperature in multi-dimensional simulations. When the 
reaction time scale becomes too long to achieve the equi- 
librium in fast flows, one needs to treat further the dy- 
namical formation or destruction of light clusters. 

Finding the light clusters in supernova cores motivates 
detailed studies on weak reactions regarding light clus- 
ters. Their reactions must be taken into account, in prin- 
ciple, in any neutrino transfer simulation of core-collapse 
supernovae and proto-neutron star cooling. It would be 
preferable to construct the supernova EOS table taking 
into account the formation of light clusters. The ap- 
proach given here merges a low-density virial expansion 
with a quasiparticle approach at higher densities. Since 
the appearance of light clusters is a natural outcome in 
dense matter, their existence has to be considered in criti- 
cal studies not only of the supernova mechanism, but also 
other problems in nuclear physics. 
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